Spectral function of the Kondo model in high magnetic fields 
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Using a recently developed perturbative renormalization group (RG) scheme, we calculate analy- 
tically the spectral function of a Kondo impurity for either large frequencies uj or large magnetic field 
B and arbitrary frequencies. For large ui 3> max[B, Tk] the spectral function decays as 1/ \v?[uj /Tk\ 
with prefactors which depend on the magnetization. The spin-resolved spectral function displays 
a pronounced peak at u ~ B with a characteristic asymmetry. In a detailed comparison with 
results from numerical renormalization group (NRG) and bare perturbation theory in next-to-leading 
logarithmic order, we show that our perturbative RG scheme is controlled by the small parameter 
1/ ln[max(o>, B)/Tk]- Furthermore, we assess the ability of the NRG to resolve structures at finite 
frequencies. 

PACS numbers: 72.15.Qm,73.63.Kv,75.20.Hr 
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I. INTRODUCTION 

Only recently it has become possible to measure 
directlyi*2iii the so-called Kondo resonance, i.e. a sharp 
resonance characteristic for the spectral function of a 
single magnetic impurity in a metal. The Kondo reso- 
nance is probably the most direct manifestation of the 
Kondo effect: The antiferromagnetic coupling between 
conduction electrons and a localized spin causes a com- 
plex many-body resonance which ultimately screens the 
magnetic degree of freedom. The spectral function is di- 
rectly proportional to a tunneling current, e.g. from the 
tip of a scanning tunneling microscope to a magnetic im- 
purity on the surface of a metallic hosli. Alternatively 
it can be measured by tunneling through weak links into 
a quantum dot in the Kondo regime^. This also opens 
the possibility of investigating the spectral function out 
of equilibrium^. 

While the spectral function is an interesting quantity 
in its own right, it also has many important applications. 
For example, it is directly related to the T-matrix, which 
describes the scattering of the conduction electrons (see 
below). It therefore determines not only the transport 
properties of magnetic impurities and quantum dots in 
the Kondo regime but it also plays an important role 
in controlling the distribution 5 of electrons in metals 
contaminated by magnetic impurities. Furthermore, the 
spectral function of impurity models is central in deter- 
mining the properties of lattice models within dynamical 
mean field theory^. 

The spectral function of the Kondo model in a 
magnetic field has been calculated using a variety of 
techniques2*LIL2ii2iiiii2iiiii4ii£, including diagrammatics, 
mean-field treatments and quantum Monte-Carlo. The 
most precise method, however, is the numerical renormal- 
ization groupie (NRG) and its generalizationsiiii^iiSiiS 
which allow the calculation of dynamical quantities. 
However, little is known analytically about the proper- 
ties of the spectral function, particularly for high mag- 
netic fields. For small B, the spectral function for small 
frequencies can be calculated from Fermi liquid theory 



and renormalized perturbation theoryi&i^. We have re- 
cently developed^ a perturbative renormalization group 
scheme based on frequency-dependent couplings for the 
Kondo model. While this method was formulated to de- 
scribe the Kondo effect out of equilibrium, it can also 
be used to calculate dynamical quantities in equilibrium 
in a controlled way for either large frequencies or large 
magnetic fields. We therefore use this method to deter- 
mine analytically the spectral function in leading order 
of 1/ In [max (B, lj)/Tk], focusing our attention on equi- 
librium and vanishing temperature. This clarifies some 
long-standing questions on the asymptotics of the spec- 
tral function in the Kondo model. A detailed compari- 
son with results from NRG and with perturbation theory 
allows us to check the assumptions underlying our per- 
turbative RG approach. 



II. THE MODEL 

We consider the spectral function ImG^(w) of an elec- 
tron in an impurity orbital within the Anderson model 

k,(T a 

+ vJ2(<iJ° + h - c -) + Un ™ n u W 
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where cj_ and /J are creation operators for electrons 
in the conduction band and in the impurity orbital, re- 
spectively, nja — flfa, and U is a large Coulomb inter- 
action matrix element. We focus our attention on the 
Kondo regime, e f < ~V 2 N F , e f + U > V 2 N F (N F 
being the conduction electron density of states), where 
the impurity orbital is occupied by a localized spin S 
(with S — 1/2) and the problem can be mapped onto the 
Kondo Hamiltonian 
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FIG. 1: Diagrammatic form of the RG equation The 
strokes symbolize derivatives with respect to In D, see Ref . I2H 
for details. 



where r is the vector of Pauli matrices. Note, that we 
have coupled an extra magnetic held B to the impurity 21 , 
measuring B in units of the Zeeman splitting. 

The Green function G*£ within the Anderson model is 
directly related to the T-matrix T a {w) = V 2 G f a {uj) of the 
conduction electrons, defined by 

sWM = AWsQl-W + AMW&W (3) 

where g c a k k , (g^) is the (bare) Green function of the con- 
duction electrons. While it is not completely obvious how 
to dehne a spectral function for a Kondo impurity, the T- 
matrix of the Kondo model can easily be identified with 
the following correlation function, e.g. using equations 
of motion^, 

T„(w) = V 2 GI(lo) = J(S Z ) + J 2 ((Sc+r aCT ; ST m ,c tt -)) , 

(4) 

where ((...)) denotes a retarded correlation function. 
Formally, we will consider the usual "scaling limit", 
where the ratios of frequency or magnetic field and 
Kondo temperature, lj/Tk, B/Tk, are kept fixed while 
the ratio of Tk and all "high energy" scales like band- 
width Do, level position or interaction are sent to zero, 
Tk I Do,Tk /U,Tk /e/ — ► 0. In this limit, the dimension- 
less product of the T-matrix and the density of states, 
iVj?ImT(w), is a universal function of uj/Tk and B/Tk 
and it coincides with the low-frequency part of the spec- 
tral function N FV 2 ImG^(u)) of the Anderson impurity. 
Therefore NfIwT^lj) is the quantity to be studied be- 
low. 



III. PERTURB ATIVE RG 

The perturbation theory of the Kondo model is charac- 
terized by logarithmic divergences. The method of choice 



for the resummation of leading logarithmic corrections is 
the perturbative renormalization group (RG). In Ref. I20I 
we have developed a certain formulation of perturbative 
RG based on frequency-dependent coupling constants. A 
thorough discussion of this approach will be published 
elsewhere 2 ^; here we will use the RG equations without 
further derivation and check the results extensively in a 
comparison with bare perturbation theory and NRG. 

The RG is formulated in terms of frequency- and spin- 
dependent vertices^ 9^ 0J '"-yJ' where a and uj c (7 and 
u)f) denote the spin and energy of the incoming elec- 
tron (the incoming spin), see Ref. l20t Primed quan- 
tities refer to outgoing particles. Before renormaliza- 

(7LJ 'Cr'uj' 

tion starts, the vertex g^J". 1 i ( J is just given by the 
bare coupling NpjT-yyiT aa i independent of frequencies. 
In leading order of the relevant expansion parameter 
1/ ln[max(w, B)/Tk], one can set the energy of the spin 
on-shell uif = —jB/2,uj^ = —j'B/2 to keep track only 
of the energy of the incoming electron. The vertex is 
parametrized as 
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)~g^- 1 B/2). (5) 



The one-loop RG equations (see Fig. ^) read2£ in terms 
of the running cut-off D 



dlnD 
dlnD 



= -2g ± (B/2) 2 Q(D - \w + <tB\) 



(6) 
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with g±(u>) = g±(—u>), g z ^{io) — g z i(—uj) and the initial 
conditions g za {iS) = g±(to) — JNp at the bare cut-off D . 
The Theta function O(w) describes that various virtual 
processes are not possible if the running cut-off D gets 
too small. Note that one recovers the usual poor man's 
scaling equations^ if one considers only the limit uo — > 
as it is usually done in RG schemes. 

The RG equations I© are easily solved and we obtain 
for D -> 



g za {w) = e[\u + <rB\ - B] 
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where within our one-loop approach Tk = Doe~ 1 ^ 2NF,J \ 
As expected, the dimensionless renormalized couplings 
are functions of ui/Tk and B/Tk only. 

It is important to realize that the logarithmic renor- 
malizations resummed in Eq. JBJ are finally cut off by 
the relaxation rate T of the spin. As the relevant pro- 
cesses involve at least one spin-flip, we identify Y with 
the transverse spin relaxation rate I/T2, which is given 
in terms of the renormalized couplings as 20 



7=-l,l 



did 



~gM 2 L(l-f w ) 



+ g ± (Lu- 1 B/2) 2 Ml-U-~ l B) ■ (8) 



where f u is the Fermi function. The rate Y cuts off 
the power-law singularities at ui = ±B (or to — ±B/2) 
and we implement this by replacing B/\tu + aB\ by 
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The errors in- 
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troduced by this heuristic replacement are small** in 
l/ln[max(.B,u;)/TK]. For B > T K , we obtain Y « 
nB/(16\ii*[B/T K ]) > T K and therefore g za (id) < 1 and 
<7_i_(u;) <C 1. Accordingly, the perturbative RG is valid 
in this regime*^ for all frequencies (while for B < Tk 
it can be used only for uj ^> Tk)- The smallness of 
the renormalized couplings for either large B or large 
uj is the main reason why our calculation is controlled by 
1/ ln[max(w, B])/Tk] as e.g. higher-order corrections to 
© remain small. 

To check the assumptions underlying our formulation 
of the perturbative RG for w-dependent vertices, we have 
calculated the vertices in next-to-leading logarithmic or- 
der in bare perturbation theory by evaluating the dia- 
grams shown in Fig. [3 Such a test seems to be useful as 
RG schemes for w-dependent vertices are not used very 
often. Neglecting terms of order g 3 InDo, we obtain 



gzo{u) 




where g — NpJ- The perturbative RG is expected to 
resum all leading logarithmic terms. Indeed, if we ex- 
pand (0 carefully in g using that Tk = -Doe -1 ' we 
obtain the same logarithmic divergences in all relevant 
limits, id — > ±B/2 or id — ► ±B, id — * ±00, and id — ► 
[the singularity at id = in (|10(l cancels]. One assump- 
tion which has not been checked in the perturbative cal- 
culation shown above is that the logarithms are cut off 
by the spin-relaxation rate. To analyze this effect one 
would have to include self-energy insertions and corre- 
sponding vertex corrections as will be discussed in a fu- 
ture publication 2 *!. 

To leading order in l/ln[max(w, B)/Tk] the spectral 
function can be calculated by replacing the bare coupling 
constants in the lowest order expression (inset of Fig. 
with the renormalized vertex (J5J and one obtains*^ 

N F lm[T a (w)] =-J(sz ff M 2 

+ 2g ± (u + a|) 2 [1 + <jM(2/ r (w + aB) - 1)]) (11) 

where M = 2(S Z ) is the magnetization of the impurity 
normalized to 1. The Fermi function fr(u>±B) has to be 
broadened by T, fr((d) = 1/2 — arctan[w/r]/7r. Diagram- 
matically this arises as the corresponding term in T(u)) 
contains a convolution with the spin-flip susceptibility 
X H (id) (see inset of Fig. |3J) and we used the fact that 
the w-dependence of the vertex g± (id) is sufficiently weak. 



For B 3> Tk one can replace M in (|llfl by 1 in leading 
order of \/\u[B /Tk\- However, to obtain the correct 
large-w behavior for arbitrary B it is essential to use the 
exact local magnetization. For B ~ Tk it is not possible 
to calculate M within our RG scheme, and we therefore 
use the magnetization obtained from the Bethe Ansatz 
(or the practically identical values from the NRG). 

Friedel's sum ruleSi relates the spectral function at id = 
exactly to the local magnetization NFlm[T a (u) = 0)] = 
-i sin 2 (So- where — = 1+ % M . As a further check of 

IT u 77 Z 

our approach, we compare this result to our perturbative 
RG. From Eqs. J7J) and (|lll) . we obtain asymptotically 
N F lm[T a {id = 0)] » -js in'ti/ifr] for B » Tk which 
is consistent with the asymptotic Bethe Ansata*^ result 
M = l-l/{2\n[B/T K ])-... 

A typical spectral function for a spin-down electron 
is shown in Fig. |3| the corresponding spin-up spectral 
function can be obtained from ImT-f(w) = ImT^(— id). 
For large frequencies the spectral function decays as 

N F Im[W} « -tt-I^^L forc^ico. (12) 
161n [td/±K\ 

A similar result, discussed below, has been obtained by 
Logan and Dickensi*^. 

For u ~ B the spectral function is characterized by 
a pronounced and highly asymmetric peak. The width 
of the left flank is determined by Y and is therefore of 
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FIG. 2: Diagrams contributing to the vertex renormalization 
up to order g 3 . Solid line: conduction electrons, dashed line: 
pseudo fermions representing the spin. See |2(il | for the corre- 
sponding Feynman rules. The (non-parquet) diagram (j) does 
not contribute to order <? 3 ln 2 Do- 
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FIG. 4: Line shape of the spectral function as a function 
of uj/B calculated within perturbative RG. All curves are 
rescaled by their maximum. For increasing B the peaks get 
more and more asymmetric. Note that a small systematic 
shift of the peak position towards lower frequencies for smaller 
B is not included (see text). Inset: peak height as a function 
olB/T K . 
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FIG. 3: T-Matrix -Npla^T^u)] (or equivalently the spec- 
tral function — NpV 2 lmG ^ (w) of an Anderson impurity) as 
a function of ui/B for B — 100 Tk- Solid line: perturbative 
RG, open circles: NRG. The NRG is not able to resolve sharp 
peaks at high frequencies. For a direct comparison, we have 
therefore convoluted (filled diamonds) the result of the per- 
turbative RG with the resolution function used in the NRG 
code. Inset: Diagram used to calculate the T-matrix Q in 
leading order using renormalized vertices. 



order Bj In 2 [B/T K ] . As T/B ~ 1/ In 2 [B/T K ] decreases 
for increasing £?, the left flank sharpens with increas- 
ing B (see Fig. . As the width of the right flank in- 
creases for large B the peak gets more and more asym- 
metric. Formally, the line shape of the right flank is 
characterized by a power-law for T <C u> — B <C B, 
see Eq. (JJJ. However, the corresponding exponents a 
of order 1/ \u[B/Tk] are so small that it is not even in 
principle possible to extract the power-law behavior— as 



a\u[B/T] ~ hx[bx(B/T K )]/ln{B/T K ) < 1. 

Within our perturbative RG, the peak in the spectral 
function is always positioned at u = B, see Fig. 0] This 
is the correct B — > oo result. A subleading shift to- 
wards lower values is not included in our scheme which 
neglects renormalization of the real part of the self- 
energy 31 . An estimate 2 ^ suggests that this effect is of 
order B/\ii[B/Tk] consistent with results from Moore 
and Weni 2 . discussed below. For B <C Tr- it has been 
shown^i^ that the peak in the spectral function is 
positioned at 2B/3. A splitting in the total spectral 
function ImTj (oj) + ImT|(w) can be observed^ for B > 
0.5 T K , where T* K w 2.8 T K is the half width at half max- 
imum (HWHM) of the Kondo resonance at T = 0. 



IV. NRG 

The results of the perturbative RG approach to the 
equilibrium T-matrix T(uj) can be compared to nu- 
merical renormalization group (NRG) results for the 
same quantity^. The NRG procedure^ consists of 
the following steps: (i) a logarithmic mesh of k-points 
k n = A~™ is introduced about the Fermi wavevector 
kp =0, and (ii) a unitary transformation of the Cka 
is performed such that = J2k c k<? is the first op- 
erator in a new basis, / n(r ,n = 0,1,..., which tridi- 
agonalizes H c = ^Z ka ^kac\ a c ka in fc-space, i.e. H c — > 

J2aJ2Zo A ~ n/2 (fl+iJn<r + h.c). The Hamiltonian © 
with the above discretized form of the kinetic energy is 
now iteratively diagonalized by defining a sequence of fi- 
nite size Hamiltonians H N = J2a,n=o A ^ n/2 (fl+iafr^ + 
h-c) + JE a ,a' S ' fla>T*'*fo,* - BS Z for > 0. For 



5 



each N, this yields the excitation energies and many- 
body eigenstates | A) jv of H N at a corresponding set of 
energy scales lun defined by the smallest scale in Hn, 

N — 1 

ljn = A 3—. Since the number of states grows as 
4 , for N > 6 only the lowest 500 or so states are re- 
tained for Hn- This limits the width of the spectrum 
of H N to < lu < K(A)u> N , where K(A) « 5 for 
the value A = 1.5 used in this paper. We are inter- 
ested in ImT(cj) = lm((0;0+)) p , where O = Sc^tw 
is defined in Eq. and the subscript p indicates the 
density matrix used to evaluate the thermodynamic av- 
erages (see below). The matrix elements (m\0\n)N, 
which are required for T(lu), are also calculated itera- 
tively. ImT(ijj) is then constructed at a characteristic set 
of frequencies uj = Sljv from Hm for each N — 0, 1, . . . 
via its Lehmann representation. This necessarily yields 
ImT(w) in the form of a set of weighted delta functions 
5{Qn ± Continuous spectra are obtained by re- 

placing the delta functions by Gaussians of width 77^: 
S(Q N ± £f ) - ^ exp(-(Q N ± E*Y/r? N ). As noted 
above, the spectrum of Hn is known in a finite win- 
dow of size ps 5uin for A = 1.5. We therefore choose 
f2jv = 1.8a; 7v, not too low since the lower part of the spec- 
trum is refined in later iterations and not too close to the 
upper bound of the spectrum which could be subject to 
truncation errors. We also set t]m — 0.39 SI jy of order the 
level structure of Hn- The main interest in this paper is 
to compare the NRG results for ImT(aj, B,T — 0) with 
the perturbative RG calculations in the limit B ^> Tk 
where the latter are expected to be accurate. In this 
limit, the spectral features in ImT(w,-B,T = 0) lie at 
high frequencies lu ~ B. A description of this finite- 
field, high-frequency part of the spectrum within NRG 
requires the use of a reduced density matrisii* 3 ^ calcu- 
lated from the true ground state, i.e. Hn-+oc, m place of 
the grand-canonical density matrix obtained from Hn- 
The necessity to use ground-state properties to calcu- 
late high-frequency features is for example evident from 
Eq. H12JI: the magnetization determines the large-o; tails. 
Nevertheless, the two approaches have been shown to 
give almost identical results at B = and to be in good 
agreement for fields B < IOTk- 



V. COMPARISONS 

How well does the perturbative RG reproduce the spec- 
tral functions calculated within NRG? Fig. shows ex- 
cellent agreement 3 for B/Tk = 100 both for lu — > and 
lu — > ±00. However, large deviations arise for lu m B 
where the peak within the perturbative RG is much 
sharper and higher. This disagreement arises because 
the logarithmic discretization in energy used within NRG 
cannot resolve sharp features at high frequencies. More 
precisely, for a feature at a given frequency Slo with in- 
trinsic width Ao it will fail to resolve i"^ if fio > Ao since 
the broadening used at frequency S7o is necessarily of or- 
der fig (0.39 Hq in our calculations 35 ). This is the case 




FIG. 5: Comparison of NRG (open symbols) and perturba- 
tive RG (filled symbols) broadened by the NRG-resolution 
function (see caption of Fig. 01 for large magnetic fields, 
5 = 1, 10,100, 1000 T K . 



here, since the spectral feature at Qq = B has a width 
on its left flank given by Ao ~ B/ In 2 [B /Tk] ■ Conse- 
quently, the width of this feature will not be captured 
by NRG for B ^> Tk- To allow a direct comparison of 
NRG and perturbative RG, we have broadened the lat- 
ter at the energy lu by convolution with the resolution 
function used in the NRG code, i.e. by a Gaussian of 
width 0.39 a;. After this broadening, the agreement of 
perturbative RG and NRG is excellent for all frequencies 
as long as B ^> Tk, see Figs. GO and Qj] For smaller B, the 
spectral function cannot be calculated reliably for small 
lu within our perturbative RG and indeed strong devia- 
tions can be seen in Fig. [S] for B < 10 Tk- Fig. shows 
that the large lu behavior is nevertheless well reproduced 
which firmly establishes Eq. I|12l) . Note that we replaced 
M in 111|) by the true ground-state magnetization for this 
comparison as discussed above. 

We have argued that the small parameter con- 
trolling the validity of our perturbative RG is 
l/ln[max(f?,a;)/Tft:]. To check this claim we plot in 
Fig. the relative error of the perturbative RG (com- 
pared to NRG) multiplied by ln[Vo; 2 + B 2 /T K ] for a 
wide range of lu and B. As expected, one obtains a 
number of order unity! Fig. [7| suggests that the next- 
order correction is suppressed by a factor smaller than 
1.5/ ln[max(i3, lu)/Tk\- From the large- B expansion of 
quantities like the magnetization calculated by Bethe 
AnsataiiSS one expects that the correction rather has the 
form ln[ln[max(_B, lu)/Tk]]/ ln[max(_B, lu)/Tk] but our 
range of parameters is too small to extract a ln[ln[. . . ]] 
prefactor. 

In the remainder of this section, we briefly compare 
our results to other approaches. Many years ago, Frota 
and Oliveiraii described the spectral function for B = 
with the simple heuristic form NfT(lu) — — 1/Re[7r(l — 
ilu/Tk ) 1//2 ] where Tk ~ 0.4 Tk (in our units) was fit- 
ted to the width of the Kondo peak. This form was 
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FIG. 6: Behavior of the spectral function — Nf-Im[Ti(ui)] = 
-N F V 2 lmGl(u)) for large lo and B = 0, WT K - For B,T K < 
lo the spectral function decays as l/ln 2 [w/Tif] with prefac- 
tors depending on the magnetization [see Eq. lilt ]. Symbols: 
NRG, solid and dot-dashed line: perturbative RG. For very 
large lo > 200 Tk one can see that an interpolation formula 
suggested in Ref. ^3 (dotted line) decaying as 1/\/lo is not 
consistent— with the NRG data. 
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lowing sense: For Tk <C lo <C D the spectral function 
decays as 1 / In 2 [lo /Tk] and not as 1 / y/u as has previ- 
ously been pointed out by Dickens and Logan 2 ** Does 
a 1/y/uJ law hold at intermediate frequencies, e.g. for 
Tk <C lo <C T*l This would require the existence of a 
new scale T* parametrically larger than Tk in contra- 
diction not only to the perturbative RG but also to the 
exact Bethe Ansatz solution 2 ^. Nevertheless, the Frota- 
Oliveira formula is a very good heuristic description in a 
large regime. For 20 Tk < lo < 200 Tk it is almost indis- 
tinguishable from our perturbative RG result for B = 
(and fits the exact result much better for lo < 20 Tk 
where our perturbative scheme breaks down). 

Recently, Logan and Dickensi^SS calculated the spec- 
tral functions in a magnetic field within their "local mo- 
ment approach" (LMA), which is based on a combination 
of diagrammatic perturbation theory and mean-field the- 
ory. Interestingly, Logan and Dickens obtain the cor- 
rect lo — ► oo and lo — > asymptotics. However the 
peak in the spectral function close to lo = B appears 
to be much broader in their approach and they find that 
the peak is shifted towards larger frequencies of order 
B h\(B /Tk)- We think that both features are an artifact 
of their heuristic scheme to introduce broadening. 

Moore and WeniS calculated the density of states of 
the spinons obtained from the Bethe Ansatz solution 28 . 
They suggested that their result may serve as an approx- 
imation to the spectral function. However, the spinon 
density of states lacks the logarithmic tails characteriz- 
ing the large lo behavior of the spectral function and their 
peaks close to lo = B do not show the correct asymmet- 
ric line shape for B 3> Tk- Nevertheless, it is interesting 
to note that the width of the spinon density of states is 
of order ~ B/ In 2 [B/Tk] like the left flank of our result. 
In their result, the position of the peak in the spectral 
function is shifted by —B/ (2 1u[B/Tk)) away from lo = B 
for B ^> Tk- As discussed above, such a shift is beyond 
the precision of our present calculation but an interesting 
open question is whether the peak-position in the spinon 
density of states coincides with that in the electron spec- 
tral function. 



FIG. 7: The perturbative RG is a controlled expansion in the 
small parameter 1/ ln[max[tj, B]/Tk]- This is shown by plot- 
ting the relative error of the perturbative RG compared to 
NRG, gg^egpg , multiplied by ln[Vw 2 + B 2 /T K ] for vari- 
ous frequencies and magnetic fields. The perturbative RG re- 
sults have been broadened by the NRG resolution (see Fig. 
and text). Only data points with \/lo 2 + B 2 > 10 Tk are 
shown. 



chosen to interpolate between the exact lo = result 
and an asymptotic l/V^ behavior which was believed 
to arise as a consequence of an X-ray singularity in the 
presence of a phase-shift 7r/2. The fit to the Frota and 
Oliveira formula works surprisingly well for frequencies 
up to 200 Tk (see Fig. |SJ| . From our point of view, this 
remarkable agreement is, however, accidental in the fol- 



VI. CONCLUSION 

In this paper we have addressed three problems in par- 
ticular. 

(i) We have calculated analytically the behavior of the 
universal spectral function in the Kondo model for either 
large frequencies or large magnetic field and arbitrary 
frequencies. We hope that our results establishes the ex- 
istence and detailed shape of the highly asymmetric peak 
at frequencies of the order of the Zeeman splitting and 
the asymptotic behavior at high frequencies. Hopefully, 
the spectral function in a magnetic field can be studied 
in some detail in future experiments e.g. by tunneling 
into magnetic impurities. 

(ii) In a detailed comparison with perturbation theory 
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and NRG, we have argued that our recently developed 
perturbative RG schemed allows controlled calculations 
in the small parameter 1/ ln[max(w, B)/Tk]- We expect 
that a similar statement holds out of equilibrium for large 
bias voltages V 3> Tk where a comparison with numeri- 
cally exact results is not yet possible. 

(iii) We have tried to assess the ability of the NRG 
to resolve structures at high frequencies. As is obvious 
from the formulation of NRG, structures which are much 
sharper than their typical frequency are not correctly re- 
produced by NRG. However, our results suggest that the 
error induced by the logarithmic discretization mainly 
results in a simple broadening. The spectral function 



for large B can serve as a test for future improvements 
of the NRG algorithms adapted to describe high-o; fea- 
tures more precisely, which may be useful, for example, 
in implementations of dynamical mean field theory using 
NRG2&. 
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